Semileptonic B a\lv decay in QCD 



T. M. Aliev * M. Savci t 
Physics Department, Middle East Technical University 
06531 Ankara, Turkey 



Abstract 

The form factors and the branching ratio of the B — > a\lv decay are calculated 
in framework of QCD sum rules. A comparison of our results on form factors and 
branching ratio with the results from constituent quark model is presented. 
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1 Introduction 



Inclusive and exclusive decays of heavy flavors play a complementary role in determina- 
tion of fundamental parameters of the Standard Model and of a deeper understanding 
the dynamics QCD. Among of all these decays the semileptonic decays occupy a special 
place, since their theoretical description is relatively simple and due to the possibility of a 
more precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
experimentally. 

The new experimental results in studying B meson decays, which have been carried 
in recent years, has improved the values of the CKM matrix elements and CP violation 
parameter. The experimental prospects at future B factories Belle [l[ and BaBar || , where 
many inclusive and exclusive channels are expected to be measured more precisely, pushes to 
make further analysis on the theoretical side. One of the main goals in these investigations 
is a more accurate determination of the CKM matrix element V u b. For this aim, both 
exclusive and inclusive b — > u transitions will be analyzed. Studying inclusive decays 
theoretically is easy, but their measurement in experiments are difficult. Using inclusive 
decays in determining V u b, implies the need of using perturbative QCD methods in the 
region near the end-point of the lepton spectrum, where many resonances are present and 
perturbative results are less reliable. This problem can be avoided by considering exclusive 
channels for which measurements in experiments are easy. But unfortunately theoretical 
description of the exclusive decays is not as easy as the inclusive decays. This is due to the 
fact that, in investigating the exclusive decays there appears the problem of calculating the 
hadronic matrix elements, which is directly related to the non-perturbative sector of QCD. 
So in an investigation of the exclusive decays, we need non-perturbative methods, such as 
QCD sum rules, lattice calculations, etc. 

The relevant semileptonic exclusive decays in determining V u b are B — > niu, B — > plv 
and B — > a\lv. The decays B — > iriv and B — >• p£v, have been extensively studied 
theoretically in a series of papers, in framework of non-perturbative approaches such as 
QCD sum rules Pj-fTI], light cone QCD sum rules [^] — (JTTJ , lattice calculations flT2|-[r3 



quark model Note that the B — > %£v and B — > plv decays have already been observed 
in experiments conducted by the CLEO Collaboration j0|. The semileptonic decay B — > 
a\lv was investigated in framework of the constituent quark model (CQM) in fll6| . Yet, this 
decay has not been observed in experiments that have been conducted so far, but expected 
to be detected in future 5-factories. For these reasons, it is the right time to investigate this 
decay in framework of different methods and a comparison of the predictions of different 
approaches is necessary. In this work we investigate the semileptonic decay B — > a\iv in 
framework of the three point QCD sum rules method. 

The present paper is organized as follows. In section 2 we formulate the sum rules for 
B — > a\ transition form factors. Section 3 is devoted to the numerical analysis of the sum 
rules, and a comparison of our results with the quark model predictions. 
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2 Sum rules for B — > a\ transition form factors 

In calculation of the B — > a\ transition form factors, we start by considering the following 
correlator 

q) = i 2 f d A xd A ye^ x -^ (o |t {d(x) lul5 u(x)J,(0)b(x) l5 d(y)}\) , (1) 



where J M = «7 M (1 — 75)6 is the weak current, with momentum q, d^y^u and bi^d are the 
currents with quantum numbers of the final a\ and initial B mesons with four-momentum 
p' and p, respectively. 

The correlator is expressed through the invariant amplitudes as follows: 



(2) 



In what follows we calculate the correlator (1) at negative values of p 2 and p' 2 with the 
help of the operator product expansion in QCD, on one side, and saturating (1) by the 
lowest meson states in the pseudoscalar and axial vector channels, on the other side. Both 
representations are equated using Borel transformations in p 2 and p' 2 , which suppress higher 
resonance and continuum contributions, as well as higher dimension operator contributions. 

Firstly let us consider the physical part of (1). Saturating (1) by the contributions 
coming from B and a\ mesons, we have 



n^(p,p') 



p 2 — m 2 B p' 2 — TOq 







d7„7 5 « Qi) (ai |u7/i(l - 75)&| B) (B bi^ b u 0) .(3) 



The weak matrix element B — > a x can be written as (q = p — p') 

(ai(e,p')l«7/x(l-75)6|5(p)) = 

2A{ ? ] - ^(m B + m ai )\/ 1 (g 2 ) + i( e 'p) (P + P% \/ 2 (g 2 ) 

™b + m 01 m B + m ai 



W) - W) 



(4) 



where m ai and £ are the mass and the four-polarization vector of the a\ meson. Note that 
form factor V 3 can be related to V\ and V 2 in the following way 



V 3 (q 2 ) = mB+mai V 1 (q 2 ) - mB ~ mai 



2m ai ' 2m ai 



V 2 (q 2 ) , 



(5) 



and with initial condition 

V 3 (0) = V o (0) . 

The vacuum-to-meson transition matrix elements are are defined in standard way, namely 



B 



bi^^d 







m% 
— Jb , 



V2m 2 ai 



9ai 



(6) 
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Secondly, let us now consider the theoretical part of the sum rules. The theoretical part of 
the sum rules is calculated by means of the operator product expansion for the correlator 
(1). Up to dimension 6, the operators are determined by the contribution of the bare 
loop, and power corrections coming from dimension-3 (ipip), dimension-4 (G 2 ) dimension- 
5 mo(V'V')) an d dimension-6 (^^) 2 operators. Our calculations show that the contributions 
coming from (G 2 ) and (ipip) 2 are negligibly small, and for this reason we don't take their 
contribution into consideration in the present analysis. 

The double dispersion relation, for the perturbative contribution, can be written as 
follows 



Tli = JJ dsds'—^ 



Pi{s,s',q 2 ) 



p 2 )(s' — p' 2 ) 



+ sub. terms , 



(7) 



where index i describes the necessary invariant amplitudes ((i — 0; +; — ; V) (see Eq. 
(2)), and pi is the corresponding spectral density. The spectral density is obtained from 
the usual Feynman integral for the bare loop by replacing 



-2rri5(p 2 ) . 



After standard calculations for the spectral densities, we have 
3rrib s'(s — 2m 2 — s' + q 2 ) 



Pv 



Po 
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4tt 2 
3ra& 

3mf, f 



A3/2 
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(ml + sq 
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:S'( S ' 



q 



1 s + 2m 2 b ) + 



2s' 
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8tt 2 ( A 3 / 2 v b/ A 5 / 2 

=F3m^(s + s' - q 2 ) =F \{s + 2m 2 b ) =F Qm 2 b s(s' + q 2 - s) =F 3sq 2 (s + s' - q 2 ) 



(8) 



where A = s + s + q 4 — 2sq — 2s' q — 2ss' is the usual triangle function. For the power 
correction contributions we get 
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+ 



2 2 
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(9) 



In Eqs. (8) and (9) subscripts within the parentheses denote the dimension of the cor- 
responding operators. Using Eqs. (3), (4), (6), (7), (8) and (9) and using double Borel 
transformation in variables p 2 and p' 2 and equating two different representations of the cor- 
relator (1) we get the following sum rules for the form factors describing B — > a\ transition: 
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(10) 



The region of integration, which appears in calculation of the perturbative contribution, is 
determined by the following inequalities: 



< s' <s' r 



o • 



m\ + 



mi 



mi 



-s < s < s , 



:iv, 



where Sq and s' Q are the continuum thresholds in the B and a\ meson channels, respectively. 
In Eq. (10), the continuum contribution is modeled as the bare loop contribution starting 
from the thresholds So and s' Q and subtracted from bare loop contribution. 

Finally, using the expressions for the form factors, we present the differential decay 
width dT/dq 2 for the B — > a\lv, which can be written in terms of the helicity amplitudes 



H±(q 2 ) 



H (q 2 ) 



(m B + m ai )Vi(q 2 ) =F 



\ l/2 (m 2 B ,m 2 ai ,q 2 ) 
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2m aiV /g^ 
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(12) 



where ±, refers to the a>i helicities. Note that the difference of the form factors V3 — Vq 
does not give any contribution to the B — > a\iv decay, since it is proportional to the lepton 
mass (in our case electron or muon). 



3 Numerical analysis 

In calculation of the form factors A(q 2 ), Vi(g 2 ), V^g 2 ) and Vo(g 2 ), we use the following 
values of the input parameters: m^ = 4.8 GeV, m ai = 1.26 GeV, g ai = 8.9 |L7|], = 
5.28 GeV, s' = 3 GeV 2 , Uip) = -(0.24 GeV) 3 pf, m 2 = (0.8 ± 0.2) GeV 2 @ (at the 



normalization point /i = 1 GeV). For the value of the leptonic decay constant, we shall 
use f B = 140 GeV. This value of f B corresponds to the case where 0(a s ) corrections are 
not taken into account (see 0, ^1|]). Since the bare loop contribution does not involve 
0(a s ) corrections, they are not taken into consideration in calculation of the leptonic decay 
constant as well, for sake of a consistent procedure. 
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The expressions for the form factors involves two independent parameters Mf and Mf . 
According to the QCD sum rules ideology, the problem is to find a region where the results 
are practically independent of Mf and Mf, while at the same time, the power corrections 
and continuum contributions remain under control. On evidence of the existing calculations 
p2f), we expect the region of stability for three-point correlator to be at values of Borel 



sec 



parameters twice as large as in the corresponding two-point functions. From an analysis of 
the two-point functions, it follows that the stability region is 4 GeV 2 < M 2 < 8 GeV 2 ((see 
20|). Therefore we find it convenient to evaluate in the range 8 GeV 2 < M 2 < 15 GeV 2 . In 



determining the range of the other independent parameter Mf , we have used the following 
relation, 

M? m 2 1 

M 2 ~ m 2 B - m\ ~ 3 ' K ' 

from which it follows that Mf ~ Mf/3. We have checked our results weak dependence on 
this ratio. 

In Figs. 1 (a)-(d), we present the Mf dependence of the form factors A(q 2 = 0), V\(q 2 = 
0), V 2 {q 2 = 0) and V (q 2 = 0) at fixed values Mf /Mf = 1/3, s' = 3 GeV 2 and at three 
different values of the threshold, s = 33 GeV 2 , s = 35 GeV 2 , s = 37 GeV 2 . From 
these figures we observe that the variations in the results, with changing sq, are about 
~ 5%. It should be noted here that, in all cases of interest, the main contribution comes 
from the dimension 3 operator {^ipj. Our final results for the values of the form factors 
at q 2 = are presented in Table 1, in which we present the CQM predictions at q 2 = as 
well. A comparison of these results yields that the magnitude of the form factor A(q 2 = 0) 
is approximately five times larger, while the magnitude of the form factors Vo(q 2 = 0) and 
Vi(q 2 = 0) are five and two times smaller, respectively, in our case, than the ones predicted 
by the CQM (the errors in form factors predicted by CQM is about 15%, see |16|]). The 



magnitude of the form factor V^g 2 = 0) is practically the same in both approaches. Figs. 
2 (a)-(d), depict the q 2 dependence of the form factors at fixed values Mf = 10 GeV 2 
and Mf = 12 GeV 2 , under the condition Mf/Mf = 1/3, and at s = 35 GeV 2 , in the 
range < q 2 < 10 GeV 2 . The reason why we consider this region of q 2 is that the non- 
perturbative contribution becomes large and hence the operator product expansion breaks 
down. In order to extend our results to the full physical region q^ ax ~ 16 GeV 2 , we 
have used the following extrapolation formula for the form factors in such a way that they 
reproduce the sum rules prediction up to 10 GeV 2 region. We have found that the best 
fits, which satisfy the above-mentioned condition, can be expressed as 

*i(0) 



F t (q 2 ) = 

1 - a F (q 2 /m 2 B ) + b F (q 2 /m 2 B ) 

The values of the constant fit parameters a F and bp, for different form factors, are listed in 
Table 1. After integrating Eq. (12) over q 2 , and using V u b = 0.032, r B = 1.56 x 10~ 12 s, one 
can obtain the decay width and branching ratio of the B — > a\iv process. These results, 
together with CQM predictions, are summarized in Table 2, which constitute the main 
results of this work. From a comparison of these results we see that, the CQM's prediction 
of branching ratio is approximately five times larger than the QCD sum rules prediction. 
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This work ap bp 


CQM pi 


A(0) 


-0.67 ±0.10 0.72 -0.20 


0.14 


Vi(0) 


-0.42 ±0.05 -0.44 0.45 


0.81 


V 2 (0) 


-0.53 ±0.05 0.45 0.13 


0.56 


V (0) 


-0.23 ±0.05 0.86 -0.38 


1.20 



Table 1: 





This work 


CQM [16] 


T + (B -> ai £u) 


5.1 x 10 5 s- 1 


(4.6 ±0.9) x 10 7 s~ l 


r_(fi -> ai £v) 


7.5 x 10 7 s- 1 


(0.98 ±0.18) x 10 8 s- 1 


T (B -> ai^i/) 


2.7 x 10 7 s- 1 


(4.0 ±0.7) x 10 8 s- 1 


-6(5 -> a^i/) 


1.6 x 10" 4 


(8.4 ± 1.6) x 10" 4 



Table 2: 
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Figure captions 



Fig. 1 The dependence of the form factors A(q 2 = 0), Vi(q 2 = 0), V 2 (q 2 = 0), and 
Vo(q 2 = 0) on Borel parameter M 2 . In all figures (a)-(d), the dotted line corresponds to 
value of threshold Sq = 33 GeV 2 , the solid line corresponds to So — 35 GeV 2 , and dash- 
dotted line corresponds to s = 37 GeV 2 , respectively. 

Fig. 2 The dependence of the form factors A(q 2 ) ,Vi(q 2 ) ,V 2 (q 2 ) and V (q 2 ) on q 2 at 
fixed value of threshold s = 35 GeV 2 . In all figures (a)-(d), the solid line corresponds to 
the fixed value of the Borel parameter M 2 = 10 GeV 2 , and the dash-dotted line corresponds 
to M 2 = 12 GeV 2 , respectively. 
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Figure 1: 
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